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A  SMALL  SCALE  THERMOSYPHON  HEAT  EXCHANGER 


Sujanani.M.,  Kiewra,  E.W.,  and  Wayner,  P.C.,  Jr. 


Department  of  Chemical  Engineering 
Rensselaer  Polytechnic  Institute 
Troy,  NY  12180-3590 


ABSTRACT 

A  small  circular  heat  transfer  cell  was  used 
to  study  the  cooling  of  a  disc  heat  source  by  a 
thin  evaporating  liquid  film.  This  small  scale 
thermosyphon  is  a  thermal  spreader  and  a  passive 
heat  sink  in  that  the  i nterf aci al ly  induced  liquid 
flow  rates  are  controlled  by  the  heat  input. 

To  determine  the  pressure  gradient,  the  evaporating 
liquid  film  profile  is  "viewed"  interferometrical ly 
through  the  objective  of  a  microscope  which  is 
equipped  with  a  scanning  microphotometer  to  measure 
its  thickness  profile.  Experimental  data  for  the 
operation  of  the  cooling  device  with  pure  fluids 
in-alkanes,  Freon  113)  are  presented.  The  combined 
use  of  an  evaporating  thin  circular  film  in  the 
shape  of  an  extended  meniscus  and  condensation 
within  a  small  volume  gives  an  effective  heat 
spreader  for  the  cooling  of  small  electronic  heat 
sources.  The  surface  heat  flux  and  heat  transfer 
coefficient  can  be  enhanced  by  increasing  the 
i nterl i ne-1 ength/area  ratio. 

NOMENCLATURE 

a  =  radius  ratio 

A^  =  area  (m?) 

A  =  modified  Hamaker  Constant  (N-m) 

B  =  London-van  der  Waals  Constant  (N-m?) 
h  =  heat  transfer  coefficient  (W/m?-K) 
h,  =  heat  of  vaporization  (kj/kg) 
kT9  =  thermal  conductivity  of  silicon  (W/m-K) 
k  =  curvature  (1/m) 

L  =  order  of  fringe 

M  =  constant  in  Equation  (4) 

me  =  evaporative  mass  flux  (kg/m?-sec) 

P  *  pressure  (N/m?) 

Pv  =  vapor  pressure  (N/m?) 

qe  =  local  heat  flux  (kW/m?) 

0  -  heat  flow  rate  (Watts) 

r  =  radius  { m ) 

£  -  resistance  to  heat  transfer  (°C/Watts) 

R  =  dimensionless  radius 

t  =  thickness  of  wafer  substrate  (m) 

T  =  temperature  (°C) 

U  =  overall  heat  transfer  coefficient  (W/m?-K) 

v  =  velocity  (m/sec) 

Greek 

6  -  film  thickness  (m) 

a  =  wavelength  of  light  (548  nm) 

n  *  refractive  index  of  liquid 

9  =  dimensionless  temperature  =  (T-T>>)/(t._j  j 


-  kinematic  viscosity  (m?'sec' 
=  viscosity  (kg/m-secl 

-  shear  stress  ( N/m?  > 


»  mass  flow  rate 
=  density  (kg/m?' 
-  surface  tension 


(kg /sec' 
( N/ml 


Subscripts 


b  =  bottom 

c  =  cavity  (rc  =  12.5  mm) 
dry  =  dry  cell 
f  =  fluid 

i  =  edge  of  heater  (r.  =  2.5  mm) 

I  'liquid 

o  =  edge  of  silicon  substrate  (r  =  38  mm) 

r  =  radial  component 

t  «  top 

vac  =  vacuum 

z  =  axial  component 

C  »  angular  component 


INTRODUCTION 


In  order  to  optimize  the  performance  of  a 
thermosyphon  and/or  a  heat  pipe,  the  following  phe¬ 
nomena  require  evaluation  and  optimization:  (i) 
fluid  flow  with  evaporation  in  ultra-thin  films 
over  a  small  surface  area  in  which  interfacial  for¬ 
ces  predominate;  (ii)  convection  of  the  resulting 
vapor;  (iii)  condensation  over  a  large  surface 
area;  and  (iv)  the  return  flow  of  the  condensate  to 
the  zone  of  evaporation.  In  this  paper  we  are  pri¬ 
marily  concerned  with  the  stability  and  the  heat 
transfer  characteristics  of  evaporating  thin  films 
which  are  two  of  the  critical  areas  to  study. 


We  present  results  concerning  the  design  and 
use  of  the  small  scale  thermosyphon  shown  in 
Figure  (1).  This  particular  design  has  both 
applied  and  experimental  attributes.  Using  the 
phenomena  associated  with  interfacial  forces  and 
the  mechanisms  of  evaporation-condensation  in  a 
heat  exchanger,  it  should  be  possible  to  transform 
the  high  heat  flux  at  the  surface  of  a  small  disc 
heat  source  to  a  lower  heat  flux  over  a  larger  area 
with  a  minimum  temperature  difference  penalty.  As 
an  experimental  tool,  the  thermosyphon  presented  in 
Figure  (1)  can  be  used  to  study  and  thereby  opti¬ 
mize  the  basic  mechanisms  of  heat,  mass,  and  momen¬ 
tum  transfer  in  an  environment  designed  to  enhance 
natural  evaporative  cooling  due  to  both  body  and 
interfacial  forces. 
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Fig.  1  Cross-Sectional  View  of  the  Small  Scale  Circular  Thermosyphon 


We  note  that  the  results  presented  below  are 
concerned  with  a  single  circular  contact  line. 
Therefore,  the  heat  fluxes  based  on  the  larger 
substrate  area  are  small.  However,  the  number  and 
length  of  interlines  on  the  substrate  can  be 
enhanced  and  therefore  the  heat  flux  can  be 
improved  by  milling,  in  a  designed  way,  multiple 
microchannels  on  the  surface.  On  the  other  hand, 
the  apparent  heat  flux  can  be  increased  by 
decreasing  the  substrate  area  relative  to  the 
interline  length.  Using  a  surface  modified  by 
capillary  grooves  and  fins,  Barthelemy,  et  al.  [1] 
have  demonstrated  high  heat  fluxes.  Herein  we  have 
reduced  the  internal  surfaces  to  a  simple  smooth 
film  which  can  be  studied  using  first  principles. 
These  "calibration  studies"  will  then  be  used  in 
future  research  to  design  enhanced  surfaces  for 
high  performance  compact  heat  exchangers. 

EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

The  small  circular  heat  cell  illustrated  in 
Figure  (i)  was  used  to  study  the  cooling  of  a  disc 
heat  source  by  a  thin  evaporating  liquid.  The 
substrate  of  the  cell  is  a  single  crystal  of 
electronic  grade  silicon  that  is  7.6  x  10*2  m  in 
diameter  and  4  x  10*4  m  thick.  The  heat  source  is 
a  resistance  heater  in  the  form  of  a  thin  platinum 
compound  film  with  an  inherent  resistance  of  about 
100  ohms.  This  S-shaped  heater  was  inscribed  in  a 
circle  of  diameter  5  x  10*3  m  on  -the  backside  of 
the  substrate,  thus  providing  intimate  contact  with 
the  silicon.  The  choice  of  using  a  single  crystal 
silicon  substrate  was  governed  by  its  high  degree 
of  surface  uniformity,  high  reflectivity  and  well- 
known  optical  and  physical  properties.  Further, 
the  processing  of  small  scale  temperature  sensors 
and  micro-grooves  planned  for  future  studies,  would 
be  possible  with  a  single-crystal  silicon  wafer  as 
the  substrate.  Uses  in  the  cooling  of  VLSI  cir¬ 
cuits  on  silicon  are  also  of  interest.  Tne 
substrate  was  epoxied  to  a  teflon  base  and  the 
outer  walls  of  the  cell  were  constructed  of 
polished  aluminum.  A  pyrex  glass  port  covered  the 
cell  and  allowed  an  unobstructed  view  of  the  heat 
transfer  chamber,  base  and  fluid.  Two  0-rings  made 
of  a  material  which  was  tested  for  inertness  with 
the  fluids  provided  a  static  seal  of  the  chamber. 

For  the  alkanes,  Buna-N  was  found  to  be  the  best 
suited  material.  The  cell  pressure  was  monitored 
and  controlled  using  a  vacuum  port.  An  injection 
port  was  also  built  into  the  cell  to  ensure  better 
control  over  fluid  volume  in  the  cell. 

To  measure  the  integral  temperature  distribu¬ 
tion  along  the  silicon  wafer,  0.0R  mm  diameter,  type 
K  thermocouples  were  attached  at  intervals  of  1  mm 


(in  the  radial  direction)  from  r  =  3  mm  to  r  =  12 
mm,  along  the  bottom  of  the  substrate.  Due  to  the 
low  thermal  conductivity  of  the  glue,  it  was  essen¬ 
tial  to  ensure  complete  and  proper  contact  of  the 
thermocouple  sensor  bob  with  the  backside  of  the 
substrate.  Since  the  thermal  conductivity  of  sili¬ 
con  is  high  (k  =  120  W/m-K) ,  the  temperature  dif¬ 
ference  across  a  thin  wafer  is  small. 

The  cell  was  assembled  and  connected  through 
the  use  of  four  clamping  screws.  Cleanliness  is 
critical  to  our  studies  as  even  a  small  amount  of 
impurities  can  radically  change  the  film  profile. 
Hence  all  the  cell  parts  were  immersed  in  absolute, 
200  proof  ethanol  and  dried  at  100°C  in  a  dust-free 
oven.  The  silicon  substrate  was  etched  with  hydro- 
flouric  acid  and  rinsed  with  ethanol,  acetone  and 
the  working  fluid,  in  that  order,  before  being 
used.  The  glass  port  was  cleaned  in  chromic  acid, 
washed  with  ethanol  and  also  dried  in  the  oven. 

The  cell  was  evacuated  by  means  of  a  vacuum 
pump  to  a  pressure  of  below  100  N/m?.  The  experi¬ 
ments  were  done  when  it  had  been  established  that 
the  vacuum  level  remained  constant  over  a  period  of 
at  least  12  hrs.  The  power  input  to  the  cell  was 
controlled  by  varying  the  voltage  across  the 
heater,  and  temperature  data  were  recorded  for 
power  inputs  of  0.5,  1,2,3,  and  4W.  To  reduce  the 
resistance  due  to  natural  convection,  a  fan  forced 
air  over  the  glass  top  during  the  experiments. 

The  working  fluid  (decane/hexane/pentane/Freon 
113)  was  then  introduced  by  opening  both  the  injec¬ 
tion  valve  and  the  vacuum  valve  to  the  cell.  The 
valve  to  the  pressure  gauge  was  closed  while  the 
fluid  was  being  injected.  When  the  valve  was 
reopened  the  pressure  rose  to  read  slightly  above 
the  vapor  pressure  of  the  fluid  at  room  tem¬ 
perature.  This  difference  was  attributed  to  the 
presence  of  a  small  amount  of  dissolved  air  in  the 
liquid.  It  is  planned  to  'degas'  the  liquid  to 
eliminate  this  dissolved  air  in  our  final  studies. 
The  volume  of  the  fluid  injected  was  controlled  so 
as  to  get  a  thin  film  within  the  operating  range  of 
power  input. 

The  temperature  profiles  were  determined  at 
the  various  power  inputs  for  the  different  working 
fluids.  These  were  compared  to  the  temperature 
profiles  obtained  when  there  was  no  fluid  in  the 

cel  1 . 

The  thickness  profile  was  measured  using  a 
scanning  microphotometer  [2].  The  fringes  observed 
due  to  the  constructive  and  destructive  inter¬ 
ference  of  light  are  as  shown  in  Figure  (2). 


«o.i  #1 
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Fig.  2  Interference  Fringes  for  Decane  at  1  W 

Wneo  these  fringes  were  scanned  using  a  photo¬ 
meter,  the  thickness  profile  '  vs.  r  can  be  calcu¬ 
lated  from 


( 2L  +  D\ 


The  solution  to  this  standard  hyperbolic 
Bessel  differential  equation  of  order  zero,  is 
given  by  [4,5] 


8flT)  = 


Cj(M)I0(Mff)  +  Ij  (M)K0(MTT) 


where  a  =  (r-j/r0) ,  and 

=  Modified  Bessel  function  of  the  first 
kind  of  order  1 k ‘ 

K|<  »  Modified  Bessel  function  of  the  second 
kind  of  order  Jk 1 

Non-dimensional  6  vs.  If  plots  were  thus 
obtained,  using  a  constant  value  of  U  =  25,  50,  and 
100  W/m^-K.  Representative  plots  of  the  experimen¬ 
tal  data  for  the  four  working  fluids  at  power 
inputs  of  0  =  1  W  and  0  -  4  W  are  shown  in  Figures 
(3)  and  (4).  We  note  that  since  the  definition  of 
U  is  arbitrary  (due  to  the  area  under  consideration 
and  definition  of  aT),  these  values  mainly  are  of 
relative  importance. 


.  L  =  0,1,2,. 


Tnis  thickness  profile  was  used  to  determine 
the  fluid  flow  and  the  heat  transfer  rate  in  the 
interline  region  using  the  momentum  balance  model. 

THERMAL  ANALYSIS 

The  analytical  solution  for  the  temperature 
profile  was  obtained  by  considering  a  circular 
silicon  plate  of  76  mm  diameter  and  0.4  mm  in 
thickness,  being  heated  in  a  central  circular 
region  of  5  mm  in  diameter,  and  cooled  by  heat  loss 
to  the  surroundings  through  different  mechanisms  on 
the  top  and  bottom  surfaces  with  a  constant  value 
of  the  overall  heat  transfer  coefficient,  II.  The 
boundary  conditions  for  the  solution  were 
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i )  T  =  Tj  at  r  =  2.5  mm  (edge  of  heater) 


li)  T  =  T  at  r  = 


i  i  i  )  -g-p  =  0  at  r  =  38  mm  (edge  of  silicon  wafer) 


Fig.  3  Temperature  vs.  Radius  for  Q  =  1W 


T-j  was  obtained  experimentally  by  extrapo¬ 
lating  the  profile  to  r  =  2.5  mm.  T  was  obtained 
experimentally  as  the  temperature  of  a  thermocouple 
in  the  air  beside  the  cell.  It  was  also  confirmed 
experimentally  that  (dT/dr)  at  the  end  of  the  sili¬ 
con  wafer  was  very  small. 

A  heat  balance  gives 

0r  =  0r+(jr  +  U.4,rdr(T-TJ  ,  or  (2) 

d/dTi  2Ur<—  —  i  . .  i 

dF  (r  d?>  =  —  (T’TJ  (3) 

To  non-d imensi onal i ze  the  problem  we  used 


#  ^ 

^  I  v 

*  4 


r  >  . 

v  * 


R"  =  —  ; 
r 

o 


;  M  =  ''2[)r0/ kt 


Using  these  definitions  in  equation  (3)  gives 


—  (IF  — )  =  m¥ 


dR  dF 


Fig.  4  Temperature  vs.  Radius  for  Q  =  4W 
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ANALYSIS  OF  RESISTANCES  to  heat  transfer 

To  optimize  the  performance  of  the  ther¬ 
mosyphon  we  need  to  evaluate  the  various  factors 
which  contribute  to  the  overall  resistance  to  heat 
transfer.  The  basic  processes  occurring  are:  (i) 
fluid  flow  with  evaporation  in  ultra-thin  films 
over  a  small  surface  area  in  which  interfacial  for¬ 
ces  predominate;  (ii)  conduction  through  the  liquid 
and  vapor;  ( i  i  i  )  convection  of  the  vapor  (diffusion 
is  important  if  non-condesi bl es  are  present); 

1 v '  condensation  over  a  large  surface  area;  (v) 
return  flow  of  the  condensate  to  the  zone  of  eva¬ 
poration;  and  (vi)  all  other  resistances  due  to 
conduction  through  the  solid  and  convection  from 
the  glass  port  to  the  surroundings. 

If  no  non-condensibles  are  present,  then  the 
resistance  due  to  (iii)  is  negligible.  We  note 
that  due  to  some  dissolved  air  in  the  liquid  and 
possible  minute  leaks  in  the  vacuum  (5  N/m^-hr)  we 
have  only  minimized  this  resistance  experimentally. 
To  completely  separate  the  various  resistances  is 
not  possible  at  present  due  to  coupling.  Data  from 
microthermistors  is  needed  for  this  purpose. 
However,  we  can  divide  the  resistances  into  two 
groups;  (i)  those  resistances  which  change  with  the 
physical  properties  of  the  fluid,  (e.g.  ,  conduction 
through  liquid,  convection  through  vapor,  fluid 
flow  and  return  flow  of  condensate);  and  (ii)  those 
resistances  which  do  not  change  when  we  change  the 
fluid,  (e.g.,  conduction  through  substrate,  glass, 
and  metal,  and  external  convection,  etc.). 

Using  a  vacuum  in  the  cell,  the  resistance 
offered  by  the  solid  part  of  the  cell  can  be  iso¬ 
lated  because  the  resistance  in  the  cavity  tends  to 
infinity.  If  we  hypothesize  that  these  two 
resistances  act  in  parallel,  we  have 


Then  by  the  approximate  values  of  ' U '  obtained 
from  Fig.  (4)  for  each  case,  we  can  get  Rf  by  using 
Equation  (R).  We  note  that  this  is  equivalent  to 
relating  the  overall  heat  transfer  coefficient  to 
the  individual  top  and  bottom  heat  transfer  coef¬ 
ficients  as 


ht  +  hb 


For  the  dry  plate  ht  -  D  and  therefore 
h. 


II 


dry  T~ 


(11) 


(12) 


If  'hb'is  assumed  constant  for  all  fluids  then 

ht  =  2(U-'.ldry)  ,  and  (13) 

„  _  1 


f  vU 


(14) 


Equation  (14)  is  exactly  equivalent  to  Eqn.  (R). 
Thus  ht  was  calculated  using  equation  (13),  and 
ht  vs.  •>„  is  ploted  in  Figure  (5)  where  Pv  is  the 
vapor  pressure  of  the  fluid  evaluated  at  T  = 

( T  i  +  TJ/?. 


(?) 


In  an  evacuated  cell,  Rf  =  - ,  and  therefore  Rvac  = 
Sc.  Now  if  Rvac  is  assumed  constant  when  Rf  is 
varied  then, 


R  R 


R  -  R 
vac 


(8) 


The  fluid  resistance  in  the  cavity  is  a  func¬ 
tion  of  various  resistances  due  to  fluid  flow,  con¬ 
densate  return,  etc.  However  we  expect  that  the 
resistance  due  to  convection  will  be  predominantly 
a  function  of  vapor  pressure.  We  note  that  the 
resistance  between  the  solid  and  surrounding  air  is 
small  due  to  the  large  sjrface  area.  Rased  on  pre¬ 
vious  theoretical  analysis  [6]  other  factors  which 
may  contribute  to  Rf  and  the  final  maximum  heat 
flux  for  one  interline  are  (ohfg/v)  and  (’5'hfg/  ). 
However  these  are  less  important  in  the  current 
studies  because  the  size  of  the  dry  spot  can  change. 

In  this  study  we  can  analyze  the  effect  of 
variation  of  vapor  pressure  on  the  overall  process 
for  a  homologous  series  of  alkanes.  Defining 


_I_ 

'JA 


(9) 


fig.  S  Heat  Transfer  Coefficient  vs.  Vapor  Pressure 

FLUID  FLOW  IN  AN  EVAPORATING  FILM  AND  THE  MOMENTUM 
BALANCE  MODEL 

From  a  basic  analysis  of  fluid  flow,  we  find 
that  the  flow  is  symmetric  with  respect  to  e  and  V.-. 

=  0.  Using  the  equation  of  continuity  and  the 
Navi er-St okes  equations  while  neglecting  vz  and 
inertia  terms,  the  following  simplified  equation 
can  be  written  for  the  r-direction 


A 


3Z  ' 


3  r 


US' 


Equation  (lb)  can  be  solved  for  the  velocity 
distribution  using  the  boundary  conditions 


z  =  0,  vr  !  0  (no  si  ip) 


(lb) 


iw 


t 


2-i  r_ 


(10) 


where  a  is  the  total  area  for  heat  transfer,  U  is 
the  overall  heat  transfer  coefficient  and  'rc‘  is 
the  radius  of  the  cavity  ( rr  =  12.5  mm). 


z  =  ■$ , 


3vr 

3Z 


Id 

dr 


(P) 


The  second  boundary  condition  has  a  signifi¬ 
cant  effect  on  the  velocity  distribution  because  of 
the  effect  of  gradients  in  temperature  and  con- 


. 


OJI 


centrat'on  (even  if  a  minute  amount  of  impurity  is 
present  1 . 

using  equations  (151  -  (17).  we  obtain 

v  =  a!!!l  r£ .  ,z]  + 1  ^  d«) 

r  „  ■  r  L  2  ‘  J  „  dr 


Tnp  normal  pressure  difference  at  the  liquid- 
iapor  interface  may  be  written  as 

O  =  p  -  -k  -  ••  .  in  nm  (IQ) 

v  .4  ' 

1*  w^r:h  curvature  may  be  approximated  as 


<  =  i  +  r,  (d 

r  rir2  V 


a“l  R  accounts  for  the  disjoining  pressure  [7], 
:or  <  lil  m  ,  J'  -  would  replace  R/  .  The  mass 


f  1  iw  ra* e  is 


The  local  value  of  the  heat  flux  is  then 

V  1  mehfg  <?e 

For  the  purpose  of  data  analysis,  we  make  the 
following  assumptions:  (i)  (dT/dr)  is  almost 
constant  in  the  interline  region  (then,  by  using 
the  pure  fluid  approximation,  (do/dr)  * 

(do/dTJ. (dT/dr)  =  constant);  and  (ii)  (di/d r)  is 
very  small  and  hence  the  curvature  may  be  repre¬ 
sented  as 


Then  by  fitting  our  experimental  6  vs.  r  pro¬ 
file  to  a  sixth  degree  polynomial  in  a  least  square 
sense  (as  in  Figure  (6)),  we  can  derive  the  mass 
flow  rate  and  the  local  heat  flux  using  Equations 
(22,25,26).  These  profiles  are  shown  in  Figure  (7). 


:»ee  polvn: 


=  V  •  2-  r  d? 

r 

0 

s’-u.;  equations  !H-?1' ,  we  get 


dk  4p  d_,, 
dr  ’  d r  '  J 


(21) 


w"p 


r  e 


1.5’ 


■1-  d-  ri T 

ir  dT  '  dr 


(22) 

(23) 


l"  onr  experiments,  'dT/dri  is  negative  and 
1-  'dr!  is  positive  for  a  pure  fluid.  Thus  the 
t emperat  ire  induced  surface  shear  stress  would 
n 1 nder  fluid  flow  towards  the  contact  line. 
However,  we  note  that  [2,6]  a  small  impurity 
a'fpcts  the  f 1 u 1 d  flow  to  a  great  extent  therefore 
if  tne  impurity  has  a  lower  vapor  pressure  and  a 
hl  ;her  surface  tension,  a  surface  tension  gradient 
d-.»>  to  a  concent  rat  i  on  gradient  can  significantly 
p-nance  flow  towards  the  contact  line.  Surface 
tension  effects  are  very  complicated. 


Tne  evaporative  mass  flux  may  now  be  obtained 
r y  using  Equation  fy*'  and  '25). 


(24! 


7  3  2 

_  dr  £R_  v  d_  .  :  dr 

r  gr  ?  dr  T"  "  dr  *  Y~  dr 


3  3  2 

?  or  ?_  dr  d_  ^  d  Jk 

1  r  -1  r  3  d  r  d  r  1“  0 


4  1,  ill  .  W  L. 

T?  ?  T  3  d  r  drdr 

dr 


Eig.  7  Mass  Flow  Rate  and  Local  Heat  Flux  vs. 

Radius 

rrom  the  fluid  flow  analysis  based  on  the 
experimental  film  thickness  vs.  radius  plot 
derived  from  thp  recorded  (ringe  pattern,  we  find 
that  the  present  model  which  does  not  take  into 
account  any  shear  stress  due  to  a  concentration 
gradient  (as  might  arise  due  to  any  impurity) 
predicts  a  very  small  amount  of  fluid  flow  away 
from  thp  contact  Imp  in  a  small  region  at  the 
beginning  of  thp  interline  and  subsequent  substan 


CONCLUSIONS 


tial  flow  towards  the  contact  line.  We  note  that 
in  our  previous  studies  we  have  found  that  flow 
towards  the  contact  line  is  always  present  and  is 
a  stable  situation.  We  therefore  do  suspect 
impurities  in  our  feed  fluid,  which  though  in  a 
very  small  concentration  (<0.0H)  could  affect  the 
profile.  Indeed  at  high  heat  fluxes,  instability 
or  breaking  up  of  the  interline  was  observed  as  is 
snown  in  Figure  (8). 


il  An  effective  thermosyphon  for  the  cooling  of 
a  disc  heat  source  was  designed. 

ii)  Integral  temperature  profiles  for  decane, 

hexane,  pentane  and  Freon  113  were  obtained. 

iiil  The  resistance  to  heat  transfer  was  found  to 
decrease  with  increasing  vapor  pressure  of 
the  fluid. 

iv)  A  preliminary  analysis  of  the  fluid  flow 
suggests  that  the  interfacial  shear  stresses 
due  to  temperature  and  concentration  gra¬ 
dients  are  critical  in  determining  the  fluid 
flow  rate. 

v)  The  heat  sink  capability  may  be  increased  by 
increasing  the  interline  length  by  making 
microchannel s. 
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Fig.  8  Instability  of  Decane  at  0  =  3W 
EXPERIMENTAL  RESULTS 

Comparison  between  the  temperature  profiles 
generated  by  the  analytical  solution  to  the 
modified  Bessel's  equation  and  the  experimental 
temperature  profiles  shows  that  the  heat  transfer 
coefficient  follows  the  trend  Unentane  > 

"freon  113  *  "hexane  >  "decane  "dry’ 

This  is  expected  because  as  the  vapor 
pressure  of  the  fluid  gets  higher,  the  resistance 
decreases  (again,  see  Figure  5). 

Average  values  for  the  overall  heat  transfer 
coefficient  and  ht  for  the  various  fluids  are 


Udry  =  25  W/m2-< 


udecane  ~  38  W/m^-K;  htdecane  "  ^ 

^hexane  =  50  W/m^-K;  hthexane  ”  W/m^-K 


Jpentane 


=  80  W/m^-K;  ht 


pentane 


=  110  W/m^-K 


uFreon  113  =  65  W/m^-K;  ht(-reon  n3  -  80  W/m^-K 

We  note  that  the  extent  of  the  interline 
region  is  only  a  fraction  of  a  millimeter.  To 
measure  the  physical  parameters  needed  to  charac¬ 
terize  and  model  heat  transfer  processes  on  a 
miniature  scale,  we  need  to  employ  devices  whose 
dimensions  are  on  the  order  of  the  process  itself. 
To  control  and  optimize  the  physical  parameters 
(e.g.,  film  thickness,  temperature  distribution, 
and  concentration)  we  need  devices  whose  sensing 
area  is  quite  small  and  sensitive.  Such  devices 
are  being  fabricated  using  miniature  solid  state 
processing  techniques,  e.g.,  photolithography, 
making  patterns  by  selective  etching. 

Since  it  is  expected  that  the  overall  heat 
sink  capability  of  the  interline  is  proportional 


a  a  i*  i  i  i  l  y  u '  l  * « c  i  'an  i  i  nr  t  nj" 

to  the  length  (,f  the  interline,  enhancing  the 
interline  length  through  surface  modification 
should  improve  the  heat  sink  capability  of  the 
thermosyphon.  We  are  considering  the  incorpora¬ 
tion  of  this  variable  in  our  future  studies. 
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